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Abstract 



We consider the model of random surfaces with extrinsic curvature term embedded into 
3d Euclidean lattice Z^. On a 3d Euclidean lattice it has equivalent representation in 
terms of transfer matrix K{Qi,Qf), which describes the propagation of loops Q. We 
study the spectrum of the transfer matrix K{Qi,Qf) on finite dimensional lattices. The 
renormalisation group technique is used to investigate phase structure of the model and 
its critical behaviour. 



1 Introduction 



Various models of random surfaces built out of triangles embedded into continuous space 
and surfaces built out of plaquettes embedded into Euclidean lattice Z''- have been 
considered in the literature fl], 0] . These models are based on area action and suffer the 
problem of non-scaling behaviour of the string tension and the dominance of branched 
polymers 0. Several studies have analyzed the physical effects produced by rigidity of the 
surface introduced by adding dimension-less extrinsic curvature term to the area action 
[^]. Comprehensive review of the work in this area up to 1997 can be found in 0. 

In this article we shall consider a model of random surfaces solely based on the with 
extrinsic curvature term embedded into 3d Euclidean lattice |Q. The corresponding 
transfer matrix has the form [^] 

K{Q^, Q2) = exp{-(3 [k{Q,) + 21{Q, A Q2) + HQ^)] }, (1) 

where Qi and Q2 are closed polygon-loops on a two-dimensional lattice, k{Q) is the 
curvature and 1{Q) is the length of the polygon-loop Q Q. This transfer matrix describes 
the propagation of the initial loop Qi to the final loop Q2. 

The spectrum of the transfer matrix which depends only on symmetric difference of 
initial and final loops Qi A Q2 

K{Qi, Q2) = exp{-2/?/(gi A Q2) }, (2) 

has been evaluated analytically in terms of correlation functions of the 2d Ising model in 
i- 

Our aim in this article is to consider the full transfer matrix (|l|) which depends not 
only on symmetric difference of loops Qi A Q2, but also on individual curvature terms 
k{Qi) + k{Q2). Because analytical solution of this problem is not known yet, we shall 
consider numerical evaluation of the transfer matrix on finite lattices in order to gain 
new insight into spectral properties of the transfer matrix in general case and critical 
properties of the model. 

In the next sections we shall numerically compute the spectrum of transfer matrix on 
finite Euclidean lattices of the sizes N -k M, where N,M = 1, 2, 3, 4. Number of loops 
on these lattices 7 = 2^*^ will be 2, 4, 16, 64, 512, 4096, 65536 and transfer matrices would 
be of the size 7x7 = 2^^^ x 2^^^ 

2 Basic Formulae 

The partition function is defined as 

Z{P)= Yl Kp{Qi,Q2)---Kp{Qr,Qi)=trK;, (3) 

{Ql,Q2,...,Qr} 

where K^^Qi, Q2) is the transfer matrix of size 7x7, defined above (|T]). Qi and Q2 are 
closed polygon-loops on a two-dimensional lattice of size N x M and 7 = 2^^^ is the 
total number of polygon- loops on a lattice T^. The transfer matrix (|l]) can be viewed 

^We shall use the word "loop" for the " polygon- loop" . 
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as describing the propagation amplitude of the polygon-loop Qi at time r to another 
polygon-loop Q2 at the time r + 1. 

The functional k{Q) is the total curvature of the polygon-loop Q which is equal to the 
number of corners of the polygon ( the vertices with self-intersection are not counted!) 
and 1{Q) is the length of Q which is equal to the number of its links. The length functional 
l{Qi A Q2) is defined as 

1{Q,) + l{Q,+i) - 2 • 1{Q, n Q,+i) = 1{Q, A Q,+i), (4) 

where the polygon-loop Qi A Q2 = Qi U Q2 \ Qi H Q2 is a union of hnks Qi U Q2 
without common links Qi fl Q2- The operation A maps two polygon-loops Qi and Q2 
into a polygon- loop Q = Qi A Q2. 

The eigenvalues of the transfer matrix K[^{Qi,Q2) define all statistical properties of 
the system and can be found as a solution of the following integral equation in the loop 
space n [|] 

^ Kf^iQ^, Q2) ^{Q2) = m ^{Qi), (5) 
{Q2} 

where "^{Q) is a function on the loop space. The eigenvalues define the partition function 
(i 

Zffj^{(3) = A5(/3, NM) + ... + A:;_i(/3, ATM), (6) 
and in the thermodynamical limit the free energy is equal to 

where Ao(/?, NM) is the largest eigenvalue defined at temperature P for the lattice of the 
size A^ X M. The correlation length is defined by the ratio of eigenvalues Ai(/5, NM) and 
Ao{(3,NM) as 

^nm{P) = -l/ln[A^iP,NM)/Ao{P,NM)], 

and grows if the eigenvalue Ai approaches to the eigenvalue Aq at some critical tempera- 
ture Pc- By the Frobenius-Perron theorem Ao(/3) is simple and Ao(/3) > Ai(/5) > A2(/3) > 

Because we can analyze numerically only finite-dimensional matrices we shall use 
renormalisation group approach to critical phenomena in order to extract critical in- 
dexes. The renormalisation group equation 

p = m (8) 

is defined through the correlation function 

kEif^ = kMM, N + M>N + M. (9) 
N + M N + M' ^ ' 

The critical point jSc is a fixed point (3* of the renormalisation group equation (|^) 

p* = m)- (10) 

The renormalisation group equation can be expanded near fixed point (3* into the 
series of the following form 
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Figure 1: The typical behaviour of the eigenvalues. These are eigenvalues of the transfer 
matrix K ClS db function of uj = e for the size 4x4. 

and the critical index u is equal to 



We shall use these formulas to define fixed point and the critical indexes. 

In the next sections we shall explicitly construct transfer matrices on finite Euclidean 
lattices of the sizes N -k M, where iV, M = 1, 2, 3, 4. Number of loops on these lattices 
^ — 2^^-'^ will be 2,4, 16,64,512,4096, 65536 and transfer matrices would be of the size 



3 General properties of the spectrum 

We shall fix the boundary spins on a lattice T^, let us say, in up direction. The convention 
is that we shall not count the boundary spins, thus if the lattice has the size 2x3 then 
it has six internal spins 0. In this setting when all spins are in up direction then in the 
dual picture we shall have just an empty loop. When one spin inside lattice is turned 
down, then in the dual picture we shall have simply one loop consisting of four bonds, a 
box-loop, and so on. 

Let us begin from the simplest lattice of the size 1x1 and then move to more com- 
plicated ones in order to understand general properties of the spectrum. For this lattice 
which has just one spin we have empty loop and one box-loop □. The transfer matri- 



ces K and K describing transitions between these two loops can be computed by using 
definition (P and (|) 



y = 



ln[{N + M) / {N + M)] 
Ingi 



(11) 



7x7 




(12) 



and then it is easy to find its eigenvalues: 



2Ao,i = 1 + e-^'^ ± ^1 + 4e-24/3 + e-ie/? - 2e-8/5, Ao,i = 1 ± 6"^^. 



(13) 



For the lattice of the size 1 x 2 we have four loops: the empty loop, two one-box loops 
^And ten boundary spins which we shall not count. 




Figure 2: The spectrum of the transfer matrices K of the size 16 x 16 as a function of 
uo = e~'^^. The quahtative character of the spectrum is the same as in Figure 1, but now 
we can see that at high temperatures some of the eigenvalues become negative. 



and one two-box loop. The transfer matrices are therefore four by four: 



/ 1 



K 



g-12/3 g-12/3 g-16/3 



e- 



g-20/3 g-16/3 



g-12/3 g-20/3 g-8/3 



-16/3 



V e 



-16/3 g-16/3 g-16/3 g- 



K 



( 1 

e" 
e" 
V e- 



e 
1 

e-12/3 1 



-12/3 



12/3 g-8/3 



-8/3 



(14) 



The eigenvalues of the matrix K have been found in 

Ao = 1 + 2e-«^ + ^ = 1 - e-^^/^, A3 = 1 - 2e-«'^ + e-^^/^, 



(15) 



and for the matrix K we have found eigenvalues on a computer (see Figure [T|). In both 
cases all eigenvalues are positive. As the last example of explicitly depicted transfer matrix 
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we shall present it for the lattice 2x2. It has sixteen loops and is of the size 16 x 16 
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(16) 



where we show only the exponents. The eigenvalues of the matrix K have been found in 



(17) 



Ao = 1 + 46-8/3 + 6e~i2/5 + 5e-i6/3, 2 3 4 = 1 + ^e-^^ - ?,e-^^^ , 
A5,6,7,8 = 1 - 2e-i2/3 + 2e-ie^ K,,^, = 1 - 2e-''P + e-^'^ 
Aii,i2,i3,i4 = 1 - 2e-8/5 + A15 = 1 - Ae-^^ + Qe-^^^ - 3e-i6/3_ 

The eigenvalues of matrix K are shown on Figure H. The general properties of the spec- 
trum at low and high temperatures can be easily understood. At /5 ^ the eigenvalues 
are : 



and when (3 



Ao = 
00 they are : 

Ao 



Ai 



A. 



2NM_l 







1, Ai 



A. 



0. 



These properties of the eigenvalues at the ends of the spectrum can be seen on Figures 1 
and 2. These Figures also allow to see the behaviour of the spectrum at the intermediary 
temperatures, the largest eigenvalue Aq and the second one Ai approach to each other. 
The rest of the eigenvalues follow the behaviour of the Ai. As we will see minimal distance 
between eigenvalue Aq and Ai decreases when the size of the lattice increases signaling 
the possibility of a phase transition. 



4 Numerical calculation of eigenvalues Aq and Ai 

To generate the transfer matrices of the larger size we have developed a special fast 
algorithm to create polygon-loops Q and then to compute the elements of the transfer 
matrix K. This part of the program is working fast enough to generate matrices of the size 
65536 and higher. To compute eigenvalues Aq and Ai we have used standard technique, 
that is acting by the matrix K on an arbitrary vector ^{Q) many, m, times. In the limit 
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m — * oo it converges to \I/o(Q)- In the same manner, acting on the orthogonal vector, 
one can get "^liQ) and therefore corresponding eigenvalues. The real problem to handle 
higher size matrices is mostly connected with the memory of the computer. 

The ratio of the eigenvalues Aq/Ai is shown on Figure |^. This ratio increases with 
the size of the matrix K and reaches the value 0.6 for the matrix of the size 4096 x 4096. 
To find fixed point j3* of the renormalization group equation we have constructed the 
curves C,nm{(3)/{N + M) for different lattices * M. As one can see on Figure 4 there 
are two fixed points at ~ 0.17 and at ^ 0.21. At the stable fixed point — 0.21 
we have computed the value gi ~ 1.99 and the critical index u ~ 0.59. This should be 
compared with the results obtained in ||13[ by low temperature expansion and by Monte- 
Carlo simulation of the corresponding system in [^, Our value is of the same order 
of magnitude. To confirm these results one should go to high-dimensional matrices. 

It is also interesting to see the behaviour of eigenf unctions. The eigenfunctions ^E'(Q) 
for the matrix K{f3) are already known and are independent of the temperature [Q. This 
is because transfer matrices K{(3) commute with each other at different temperatures and 
the system is integrable [|]. This is a general property of integrable systems |^, [1^, 0, |T^. 
The situation with the full transfer matrices K{(3) is not known, therefore we should check 
if the eigenfunctions are temperature dependent. For that reason we shall compute the 
scalar product 

between eigenfunction \E'a(<5)(/3) at different temperatures P' and /? in order to see if 
the scalar product depends on temperature. For the first eigenvalue Aq and /?' = the 
result of the calculations is shown on Figure 5. It is clearly seen that the scalar product 
is temperature dependent and that it drastically changes near the critical points. This 
result drives us to the conclusion that full transfer matrix represents a non-integrable 
system. 
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Figure 5: Temperature dependents of the scalar product {\&o(Q)(0) *^o(<5)(/5)} for the 
matrices of the sizes 16 x 16, 64 x 16, 512 x 512, 4096 x 4096. 
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